Epigenetic modulators link mitochondrial redox homeostasis to cardiac function in a sex-dependent manner.

While excessive production of reactive oxygen species (ROS) is a characteristic hallmark of numerous diseases, clinical approaches that ameliorate oxidative stress have been unsuccessful. Here, utilizing multi-omics, we demonstrate that in cardiomyocytes, mitochondrial isocitrate dehydrogenase (IDH2) constitutes a major antioxidative defense mechanism. Paradoxically reduced expression of IDH2 associated with ventricular eccentric hypertrophy is counterbalanced by an increase in the enzyme activity. We unveil redox-dependent sex dimorphism, and extensive mutual regulation of the antioxidative activities of IDH2 and NRF2 by a feedforward network that involves 2-oxoglutarate and L-2-hydroxyglutarate and mediated in part through unconventional hydroxy-methylation of cytosine residues present in introns. Consequently, conditional targeting of ROS in a murine model of heart failure improves cardiac function in sex- and phenotype-dependent manners. Together, these insights may explain why previous attempts to treat heart failure with antioxidants have been unsuccessful and open new approaches to personalizing and, thereby, improving such treatment.

Identification of microbial pathogens in Neolithic Scandinavian humans.

With the Neolithic transition, human lifestyle shifted from hunting and gathering to farming. This change altered subsistence patterns, cultural expression, and population structures as shown by the archaeological/zooarchaeological record, as well as by stable isotope and ancient DNA data. Here, we used metagenomic data to analyse if the transitions also impacted the microbiome composition in 25 Mesolithic and Neolithic hunter-gatherers and 13 Neolithic farmers from several Scandinavian Stone Age cultural contexts. Salmonella enterica, a bacterium that may have been the cause of death for the infected individuals, was found in two Neolithic samples from Battle Axe culture contexts. Several species of the bacterial genus Yersinia were found in Neolithic individuals from Funnel Beaker culture contexts as well as from later Neolithic context. Transmission of e.g. Y. enterocolitica may have been facilitated by the denser populations in agricultural contexts.

aMeta: an accurate and memory-efficient ancient metagenomic profiling workflow.

Analysis of microbial data from archaeological samples is a growing field with great potential for understanding ancient environments, lifestyles, and diseases. However, high error rates have been a challenge in ancient metagenomics, and the availability of computational frameworks that meet the demands of the field is limited. Here, we propose aMeta, an accurate metagenomic profiling workflow for ancient DNA designed to minimize the amount of false discoveries and computer memory requirements. Using simulated data, we benchmark aMeta against a current state-of-the-art workflow and demonstrate its superiority in microbial detection and authentication, as well as substantially lower usage of computer memory.

The microbial biodiversity at the archeological site of Tel Megiddo (Israel).

Introduction: The ancient city of Tel Megiddo in the Jezreel Valley (Israel), which lasted from the Neolithic to the Iron Age, has been continuously excavated since 1903 and is now recognized as a World Heritage Site. The site features multiple ruins in various areas, including temples and stables, alongside modern constructions, and public access is allowed in designated areas. The site has been studied extensively since the last century; however, its microbiome has never been studied. We carried out the first survey of the microbiomes in Tel Megiddo. Our objectives were to study (i) the unique microbial community structure of the site, (ii) the variation in the microbial communities across areas, (iii) the similarity of the microbiomes to urban and archeological microbes, (iv) the presence and abundance of potential bio-corroding microbes, and (v) the presence and abundance of potentially pathogenic microbes. Methods: We collected 40 swab samples from ten major areas and identified microbial taxa using next-generation sequencing of microbial genomes. These genomes were annotated and classified taxonomically and pathogenetically. Results: We found that eight phyla, six of which exist in all ten areas, dominated the site (>99%). The relative sequence abundance of taxa varied between the ruins and the sampled materials and was assessed using all metagenomic reads mapping to a respective taxon. The site hosted unique taxa characteristic of the built environment and exhibited high similarity to the microbiome of other monuments. We identified acid-producing bacteria that may pose a risk to the site through biocorrosion and staining and thus pose a danger to the site's preservation. Differences in the microbiomes of the publicly accessible or inaccessible areas were insignificant; however, pathogens were more abundant in the former. Discussion: We found that Tel Megiddo combines microbiomes of arid regions and monuments with human pathogens. The findings shed light on the microbial community structures and have relevance for bio-conservation efforts and visitor health.

Historical RNA expression profiles from the extinct Tasmanian tiger.

Paleogenomics continues to yield valuable insights into the evolution, population dynamics, and ecology of our ancestors and other extinct species. However, DNA sequencing cannot reveal tissue-specific gene expression, cellular identity, or gene regulation, which are only attainable at the transcriptional level. Pioneering studies have shown that useful RNA can be extracted from ancient specimens preserved in permafrost and historical skins from extant canids, but no attempts have been made so far on extinct species. We extract, sequence, and analyze historical RNA from muscle and skin tissue of a ∼130-year-old Tasmanian tiger (Thylacinus cynocephalus) preserved in desiccation at room temperature in a museum collection. The transcriptional profiles closely resemble those of extant species, revealing specific anatomical features such as slow muscle fibers or blood infiltration. Metatranscriptomic analysis, RNA damage, tissue-specific RNA profiles, and expression hotspots genome-wide further confirm the thylacine origin of the sequences. RNA sequences are used to improve protein-coding and noncoding annotations, evidencing missing exonic loci and the location of ribosomal RNA genes while increasing the number of annotated thylacine microRNAs from 62 to 325. We discover a thylacine-specific microRNA isoform that could not have been confirmed without RNA evidence. Finally, we detect traces of RNA viruses, suggesting the possibility of profiling viral evolution. Our results represent the first successful attempt to obtain transcriptional profiles from an extinct animal species, providing thought-to-be-lost information on gene expression dynamics. These findings hold promising implications for the study of RNA molecules across the vast collections of natural history museums and from well-preserved permafrost remains.

Evolutionary consequences of genomic deletions and insertions in the woolly mammoth genome.

Woolly mammoths had a set of adaptations that enabled them to thrive in the Arctic environment. Many mammoth-specific single nucleotide polymorphisms (SNPs) responsible for unique mammoth traits have been previously identified from ancient genomes. However, a multitude of other genetic variants likely contributed to woolly mammoth evolution. In this study, we sequenced two woolly mammoth genomes and combined these with previously sequenced mammoth and elephant genomes to conduct a survey of mammoth-specific deletions and indels. We find that deletions are highly enriched in non-coding regions, suggesting selection against structural variants that affect protein sequences. Nonetheless, at least 87 woolly mammoth genes contain deletions or indels that modify the coding sequence, including genes involved in skeletal morphology and hair growth. These results suggest that deletions and indels contributed to the unique phenotypic adaptations of the woolly mammoth, and were potentially critical to surviving in its natural environment.

High-throughput muscle fiber typing from RNA sequencing data.

BACKGROUND: Skeletal muscle fiber type distribution has implications for human health, muscle function, and performance. This knowledge has been gathered using labor-intensive and costly methodology that limited these studies. Here, we present a method based on muscle tissue RNA sequencing data (totRNAseq) to estimate the distribution of skeletal muscle fiber types from frozen human samples, allowing for a larger number of individuals to be tested. METHODS: By using single-nuclei RNA sequencing (snRNAseq) data as a reference, cluster expression signatures were produced by averaging gene expression of cluster gene markers and then applying these to totRNAseq data and inferring muscle fiber nuclei type via linear matrix decomposition. This estimate was then compared with fiber type distribution measured by ATPase staining or myosin heavy chain protein isoform distribution of 62 muscle samples in two independent cohorts (n = 39 and 22). RESULTS: The correlation between the sequencing-based method and the other two were rATPas = 0.44 [0.13-0.67], [95% CI], and rmyosin = 0.83 [0.61-0.93], with p = 5.70 × 10-3 and 2.00 × 10-6, respectively. The deconvolution inference of fiber type composition was accurate even for very low totRNAseq sequencing depths, i.e., down to an average of ~ 10,000 paired-end reads. CONCLUSIONS: This new method ( https://github.com/OlaHanssonLab/PredictFiberType ) consequently allows for measurement of fiber type distribution of a larger number of samples using totRNAseq in a cost and labor-efficient way. It is now feasible to study the association between fiber type distribution and e.g. health outcomes in large well-powered studies.

App-based COVID-19 syndromic surveillance and prediction of hospital admissions in COVID Symptom Study Sweden.

The app-based COVID Symptom Study was launched in Sweden in April 2020 to contribute to real-time COVID-19 surveillance. We enrolled 143,531 study participants (≥18 years) who contributed 10.6 million daily symptom reports between April 29, 2020 and February 10, 2021. Here, we include data from 19,161 self-reported PCR tests to create a symptom-based model to estimate the individual probability of symptomatic COVID-19, with an AUC of 0.78 (95% CI 0.74-0.83) in an external dataset. These individual probabilities are employed to estimate daily regional COVID-19 prevalence, which are in turn used together with current hospital data to predict next week COVID-19 hospital admissions. We show that this hospital prediction model demonstrates a lower median absolute percentage error (MdAPE: 25.9%) across the five most populated regions in Sweden during the first pandemic wave than a model based on case notifications (MdAPE: 30.3%). During the second wave, the error rates are similar. When we apply the same model to an English dataset, not including local COVID-19 test data, we observe MdAPEs of 22.3% and 19.0% during the first and second pandemic waves, respectively, highlighting the transferability of the prediction model.

Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D.

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

Increased burden of ultra-rare structural variants localizing to boundaries of topologically associated domains in schizophrenia.

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

Plasma Lipidome and Prediction of Type 2 Diabetes in the Population-Based Malmö Diet and Cancer Cohort.

OBJECTIVE: Type 2 diabetes mellitus (T2DM) is associated with dyslipidemia, but the detailed alterations in lipid species preceding the disease are largely unknown. We aimed to identify plasma lipids associated with development of T2DM and investigate their associations with lifestyle. RESEARCH DESIGN AND METHODS: At baseline, 178 lipids were measured by mass spectrometry in 3,668 participants without diabetes from the Malmö Diet and Cancer Study. The population was randomly split into discovery (n = 1,868, including 257 incident cases) and replication (n = 1,800, including 249 incident cases) sets. We used orthogonal projections to latent structures discriminant analyses, extracted a predictive component for T2DM incidence (lipid-PCDM), and assessed its association with T2DM incidence using Cox regression and lifestyle factors using general linear models. RESULTS: A T2DM-predictive lipid-PCDM derived from the discovery set was independently associated with T2DM incidence in the replication set, with hazard ratio (HR) among subjects in the fifth versus first quintile of lipid-PCDM of 3.7 (95% CI 2.2-6.5). In comparison, the HR of T2DM among obese versus normal weight subjects was 1.8 (95% CI 1.2-2.6). Clinical lipids did not improve T2DM risk prediction, but adding the lipid-PCDM to all conventional T2DM risk factors increased the area under the receiver operating characteristics curve by 3%. The lipid-PCDM was also associated with a dietary risk score for T2DM incidence and lower level of physical activity. CONCLUSIONS: A lifestyle-related lipidomic profile strongly predicts T2DM development beyond current risk factors. Further studies are warranted to test if lifestyle interventions modifying this lipidomic profile can prevent T2DM.

Spatially and functionally distinct subclasses of breast cancer-associated fibroblasts revealed by single cell RNA sequencing.

Cancer-associated fibroblasts (CAFs) are a major constituent of the tumor microenvironment, although their origin and roles in shaping disease initiation, progression and treatment response remain unclear due to significant heterogeneity. Here, following a negative selection strategy combined with single-cell RNA sequencing of 768 transcriptomes of mesenchymal cells from a genetically engineered mouse model of breast cancer, we define three distinct subpopulations of CAFs. Validation at the transcriptional and protein level in several experimental models of cancer and human tumors reveal spatial separation of the CAF subclasses attributable to different origins, including the peri-vascular niche, the mammary fat pad and the transformed epithelium. Gene profiles for each CAF subtype correlate to distinctive functional programs and hold independent prognostic capability in clinical cohorts by association to metastatic disease. In conclusion, the improved resolution of the widely defined CAF population opens the possibility for biomarker-driven development of drugs for precision targeting of CAFs.

N1-methylnicotinamide is a signalling molecule produced in skeletal muscle coordinating energy metabolism.

Obesity is a major health problem, and although caloric restriction and exercise are successful strategies to lose adipose tissue in obese individuals, a simultaneous decrease in skeletal muscle mass, negatively effects metabolism and muscle function. To deeper understand molecular events occurring in muscle during weight-loss, we measured the expressional change in human skeletal muscle following a combination of severe caloric restriction and exercise over 4 days in 15 Swedish men. Key metabolic genes were regulated after the intervention, indicating a shift from carbohydrate to fat metabolism. Nicotinamide N-methyltransferase (NNMT) was the most consistently upregulated gene following the energy-deficit exercise. Circulating levels of N1-methylnicotinamide (MNA), the product of NNMT activity, were doubled after the intervention. The fasting-fed state was an important determinant of plasma MNA levels, peaking at ~18 h of fasting and being lowest ~3 h after a meal. In culture, MNA was secreted by isolated human myotubes and stimulated lipolysis directly, with no effect on glucagon or insulin secretion. We propose that MNA is a novel myokine that enhances the utilization of energy stores in response to low muscle energy availability. Future research should focus on applying MNA as a biomarker to identify individuals with metabolic disturbances at an early stage.

Glucose-Induced Changes in Gene Expression in Human Pancreatic Islets: Causes or Consequences of Chronic Hyperglycemia.

Dysregulation of gene expression in islets from patients with type 2 diabetes (T2D) might be causally involved in the development of hyperglycemia, or it could develop as a consequence of hyperglycemia (i.e., glucotoxicity). To separate the genes that could be causally involved in pathogenesis from those likely to be secondary to hyperglycemia, we exposed islets from human donors to normal or high glucose concentrations for 24 h and analyzed gene expression. We compared these findings with gene expression in islets from donors with normal glucose tolerance and hyperglycemia (including T2D). The genes whose expression changed in the same direction after short-term glucose exposure, as in T2D, were considered most likely to be a consequence of hyperglycemia. Genes whose expression changed in hyperglycemia but not after short-term glucose exposure, particularly those that also correlated with insulin secretion, were considered the strongest candidates for causal involvement in T2D. For example, ERO1LB, DOCK10, IGSF11, and PRR14L were downregulated in donors with hyperglycemia and correlated positively with insulin secretion, suggesting a protective role, whereas TMEM132C was upregulated in hyperglycemia and correlated negatively with insulin secretion, suggesting a potential pathogenic role. This study provides a catalog of gene expression changes in human pancreatic islets after exposure to glucose.

LoFtool: a gene intolerance score based on loss-of-function variants in 60 706 individuals.

Motivation: Depletion of loss-of-function (LoF) mutations may provide a rank of genic functional intolerance and consequently susceptibility to disease. Results: Here we have studied LoF mutations in 60 706 unrelated individuals and show that the most intolerant quartile of ranked genes is enriched in rare and early onset diseases and explains 87% of de novo haploinsufficient OMIM mutations, 17% more than any other gene scoring tool. We detected particular enrichment in expression of the depleted LoF genes in brain (odds ratio = 1.5; P -value = 4.2e-07). By searching for de novo haploinsufficient mutations putatively associated with neurodevelopmental disorders in four recent studies, we were able to explain 81% of them. Taken together, this study provides a novel gene intolerance ranking system, called LoFtool, which may help in ranking genes of interest based on their LoF intolerance and tissue expression. Availability and implementation: The LoFtool gene scores are available in the Supplementary data . Contact: joaofadista@gmail.com. Supplementary information: Supplementary data are available at Bioinformatics online.

The Genetic Landscape of Renal Complications in Type 1 Diabetes.

Diabetes is the leading cause of ESRD. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2843 subjects, we estimated that the heritability of diabetic kidney disease was 35% (P=6.4×10-3). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole-exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index (P=2.2×10-5) and the risk of type 2 diabetes (P=6.1×10-4) associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation (P=1.1×10-4). Pathway analysis implicated ascorbate and aldarate metabolism (P=9.0×10-6), and pentose and glucuronate interconversions (P=3.0×10-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.

Large-scale GWAS identifies multiple loci for hand grip strength providing biological insights into muscular fitness.

Hand grip strength is a widely used proxy of muscular fitness, a marker of frailty, and predictor of a range of morbidities and all-cause mortality. To investigate the genetic determinants of variation in grip strength, we perform a large-scale genetic discovery analysis in a combined sample of 195,180 individuals and identify 16 loci associated with grip strength (P<5 × 10-8) in combined analyses. A number of these loci contain genes implicated in structure and function of skeletal muscle fibres (ACTG1), neuronal maintenance and signal transduction (PEX14, TGFA, SYT1), or monogenic syndromes with involvement of psychomotor impairment (PEX14, LRPPRC and KANSL1). Mendelian randomization analyses are consistent with a causal effect of higher genetically predicted grip strength on lower fracture risk. In conclusion, our findings provide new biological insight into the mechanistic underpinnings of grip strength and the causal role of muscular strength in age-related morbidities and mortality.

Identification of proliferative and mature ß-cells in the islets of Langerhans.

Insulin-dependent diabetes is a complex multifactorial disorder characterized by loss or dysfunction of ß-cells. Pancreatic ß-cells differ in size, glucose responsiveness, insulin secretion and precursor cell potential; understanding the mechanisms that underlie this functional heterogeneity might make it possible to develop new regenerative approaches. Here we show that Fltp (also known as Flattop and Cfap126), a Wnt/planar cell polarity (PCP) effector and reporter gene acts as a marker gene that subdivides endocrine cells into two subpopulations and distinguishes proliferation-competent from mature ß-cells with distinct molecular, physiological and ultrastructural features. Genetic lineage tracing revealed that endocrine subpopulations from Fltp-negative and -positive lineages react differently to physiological and pathological changes. The expression of Fltp increases when endocrine cells cluster together to form polarized and mature 3D islet mini-organs. We show that 3D architecture and Wnt/PCP ligands are sufficient to trigger ß-cell maturation. By contrast, the Wnt/PCP effector Fltp is not necessary for ß-cell development, proliferation or maturation. We conclude that 3D architecture and Wnt/PCP signalling underlie functional ß-cell heterogeneity and induce ß-cell maturation. The identification of Fltp as a marker for endocrine subpopulations sheds light on the molecular underpinnings of islet cell heterogeneity and plasticity and might enable targeting of endocrine subpopulations for the regeneration of functional ß-cell mass in diabetic patients.

Excess maternal transmission of variants in the THADA gene to offspring with type 2 diabetes.

AIMS/HYPOTHESIS: Genome-wide association studies (GWAS) have identified more than 65 genetic loci associated with risk of type 2 diabetes. However, the contribution of distorted parental transmission of alleles to risk of type 2 diabetes has been mostly unexplored. Our goal was therefore to search for parent-of-origin effects (POE) among type 2 diabetes loci in families. METHODS: Families from the Botnia study (n = 4,211, 1,083 families) were genotyped for 72 single-nucleotide polymorphisms (SNPs) associated with type 2 diabetes and assessed for POE on type 2 diabetes. The family-based Hungarian Transdanubian Biobank (HTB) (n = 1,463, >135 families) was used to replicate SNPs showing POE. Association of type 2 diabetes loci within families was also tested. RESULTS: Three loci showed nominal POE, including the previously reported variants in KCNQ1, for type 2 diabetes in families from Botnia (rs2237895: p POE = 0.037), which can be considered positive controls. The strongest POE was seen for rs7578597 SNP in the THADA gene, showing excess transmission of the maternal risk allele T to diabetic offspring (Botnia: p POE = 0.01; HTB p POE = 0.045). These data are consistent with previous evidence of allelic imbalance for expression in islets, suggesting that the THADA gene can be imprinted in a POE-specific fashion. Five CpG sites, including those flanking rs7578597, showed differential methylation between diabetic and non-diabetic donor islets. CONCLUSIONS/INTERPRETATION: Taken together, the data emphasise the need for genetic studies to consider from which parent an offspring has inherited a susceptibility allele.